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SECTION  I 
INTRODUCTION 

The  usefulness  of  any  constitutive  equation  depends  on  three 
major  factors:  (1)  how  well  the  constitutive  equation  can  predict 
the  mechanical  response  of  material  to  a  broad  range  of  load 
histories;  (2)  how  easily  the  equation  can  be  used  in  computer 
algorithms  and,  in  particular,  its  adaptability  to  finite  element 
analyses  of  structural  components;  and  (3)  how  easily  the  paramete 
in  the  constitutive  model  can  be  determined  from  experimental 
data.  Any  equation  that  does  not  satisfy  all  the  above  conditions 
is  severely  limited  for  its  scope  of  application. 

This  paper  is  concerned  with  procedures  for  the  dot cnu i nat  i on 
of  evolution  equations  in  a  state  variable  constitutive  thetny, 
and  a  partial  evaluation  of  the  predictive  capability  of  those, 
equations  with  the  determined  constants.  The  equations  arc 
characterized  by  two  internal  state  variables  which  roiiia'scnt  tin 
resistance  of  the  material  to  plastic  flow  and  internal  damace. 

The  particular  set  of  equations  used  in  this  exercise  an'  those 
developed  by  Bodner  and  Partom,  Refs.  [1,2],  but  other  equations 
of  the  same  general  class  could  be  treated  by  similar  methods. 

The  method  of  developing  the  state  variable  evolution  equal lo 
suggested  in  this  paper  permits  direct  evaluation  of  the  material 
parameters  from  the  experimental  data  and  thereby  minimizes 
subsequent  trial  and  error  computation.  Since  state  variable 
constitutive  equations  predict  the  strain  rate  as  a  function 
(not  functional)  of  the  stress,  temperature  and  state  variable 
histories;  the  constitutive  function  may  be  algebraically 
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inverted  to  calculate  the  state  variable  .i^story  from  stress,  strain 
rate  and  temperature  histories  that  are  determined  from  an  experiment. 

A  systematic  analysis  of  the  experimental  data  can  therefore  contribute 
to  the  formulation  of  the  evolution  equations  for  the  state  variables 
as  well  as  defining  the  material  constants. 

An  interesting  consequence  of  this  technique  is  that  it  can  be 
used  to  evaluate  more  than  one  state  variable.  This  allowed  for  the 
consideration  of  a  damage  state  variable.  Ref.  [3].  Although  a  general 
development  for  the  damage  evolution  equation  is  not  included;  a  few 
properties  of  the  variable  were  deduced  from  the  above  analysis. 

The  procedure  is  demonstrated  for  wrought  Rene  95  at  650°C  using 
data  currently  available  in  the  literature.  Refs.  [4,  5,  6].  The 
results  are  also  compared  to  am  earlier  application  of  the  Bodner  Partom 
equation.  Ref.  [7],  for  the  same  material.  Examples  of  the  calculated 
response  are  compared  to  observed  tensile  and  creep  results. 
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SECTION  II 
THE  FLOW  EQUATION 


The  basic  equations  of  Bodner-Partom  will  be  summarized 
for  reference.  These  equations  were  motivated  by  the  concepts 
and  equations  of  "dislocation  dynamics",  e.g..  Refs.  [8,  9], 
and  were  formulated  in  the  context  of  continuum  mechanics , 

Ref.  [10].  For  small  strains,  the  strain  rates  are  considered 
to  be  decomposable  into  elastic  (reversible)  and  inelastic 
(non-reversible)  components. 


e  .  . 
1] 


•  e  -p 

e  .  .  +  e  .  . 
3-1  3-D 


(1) 


*  ^ 

where  are  given  in  terms  of  stress  rates  determined  from 
the  time  derivative  of  Hooke's  Law.  Both  components  in  Eq.  (1) 
are  always  non-zero  for  non-zero  stresses  and  stress  rates 
which  imply  that  a  yield  criteria  and  loading/unloading  conditions 
are  not  required. 

The  plastic  strain  rate  is  written  in  a  form  similar  to 
the  Prandtl-Reuse  flow  law. 


(2) 


where  e? .  and  s. .  are  the  deviatoric  strain  rate  and  stress 
3-D  3.D 

tensors,  respectively.  Equation  (2)  satisfies  the  condition  in¬ 
compressibility.  Squaring  Eq.  (2)  leads  to 

dP  =  J2  (3) 

where  and  ajre  the  second  invariants  of  the  deviatoric  plastic 
strain  rate  and  stress  tensors. 
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A  basic  assumption  of  the  formulation  of  Bodner-Partom  is 
that  all  inelastic  deformations  are  controlled  by  a  relationshic 
between  d|  and  which  involves  other  state  variables,  i.e. 

d|  =  f{J^,  T,  Z.)  (4) 


where  T  is  the  temperature  and  are  history  dependent  internal 
state  variables.  Combining  Eqs.  (2),  (3)  and  (4)  leads  to  a 
representation  for  the  inelastic  strain  rate  tensor  that  can  be 
written  as  a  function  of  stress,  temperature  and  internal  state 
variables;  i.e. 


(5) 


A  particular  form  for  Eq.  (4)  was  established  from  the 
relationship  between  dislocation  velocity  and  stress,  Refs.  [8,  9], 
by  Bodner  and  Partom.  In  a  constant  temperature  environment 
this  representation  is 

f^^^  =  exp[-(Z^/3J2)”  {n+l)/(n)]  .  (6) 


In  Eq.  (6),  is  the  limiting  strain  rate  in  shear,  n  is  a 
material  constant  that  controls  strain  rate  sensitivity  and  also 
influences  the  overall  level  of  the  stress-strain  curves,  and  Z  is 
interpreted  as  an  internal  state  variable  that  governs  resistance 
to  plastic  flow.  That  is,  an  increase  of  Z  would  correspond  to 
work  hardening  and  would  require  an  increase  of  stress  to  maintain 
a  given  plastic  strain  rate.  Physically,  Z  can  be  interpretated 
as  a  measure  of  the  stored  energy  of  cold  work.  Refs.  [11/  12] . 
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In  the  earlier  references  it  was  proposed  that  Z  should  be 
taken  to  be  a  scalar  function  of  plastic  work  for  isotropic  hardening 
conditions.  This  was  later  extended  to  include  both  isotropic  and 
kinematic  hardening  by  using  different  values  of  Z  for  positive  and 
negative  stress  in  uniaxial  reversed  loading  conditions,  Ref.  [12]. 

The  model  was  subsequently  generalized  by  Stouffer  and  Bodner,  Ref. 
[13] ,  to  deformation  induced  anisotropic  flow.  In  this  case  Z  was 
shown  to  be  a  second  order  tensor  in  a  six  dimensional  vector  space. 
The  discussion  in  this  paper  is  limited  to  isotropic  hardening  for 
which  Z  is  a  single  valued  scalar  for  reversed  plastic  flow. 

Let  us  also  briefly  note  some  recent  work,  Ref.  [3],  where 
isotropic  damage  was  included  as  an  additional  state  variable  in  the 
flow  equations.  Using  the  Kachanov  definition.  Ref.  [14],  let  ui  be 
a  history  dependent  variable  that  represents  the  deterioration  in  the 
ability  of  a  material  to  support  a  stress.  In  this  formulation  stress 
j ,  is  replaced  by  a  pseudo  stress  a^^/(l-x).  Introducing  this 
result  into  the  flow  equations  ultimately  reduces  to  replacing  Z  by 
Z  (l-co). 

Applying  the  preceding  model  to  the  case  of  an  uniaxial  stress 
a,  the  one  dimensional  flow  equation  can  be  written  as 

eP  =  ^  (^)  D  exp  f-  j 

/3  1 0 1  °  2  c  n 

for  the  uniaxial  inelastic  strain  rate,  A  direct  method  of 

developing  the  evolution  equations  for  the  state  variables  Z  and 
j  is  introduced  in  the  next  section.  The  method  is  general  and  can 
be  applied  to  any  state  variable  equation.  Equation  (7)  will  be  used 
in  this  particular  example. 
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SECTION  III 


STRUCTURE  OF  THE  EVOLUTION  EQUATIONS 


Let  us  consider  the  evolution  equation  for  the  state 
variable  Z  in  the  absence  of  damage.  Therefore  tertiary  creep 
and  the  large  strain  range  of  a  tensile  test  are  excluded  from 
this  analysis.  The  evolution  equation  for  the  state  variable  Z 
is  generally  sought  in  the  form  of  a  differential  equation; 

Z  =  Z)  .  (8) 


A  more  specific  representation  is  based  on  the  concept  that 

plastic  work,  W^,  controls  the  hardening  process,  and  that 

and  are  functions  of  o.  .  and  e .  .  which  can  be  determined 

ID 

from  and  Z.  That  is,  assume 


Z  *  fj^  (W^) 


(9) 


where  f^  and  f2  are  material  functions  and  (  ) '  denotes  differ¬ 
entiation  with  respect  to  plastic  work.  Equation  (9)  is  a 
generalization  of  the  specific  representation  developed  by  Bodner. 

The  underlying  philosophy  in  assuming  the  structure  of  Eq.  (9)  is  to 
allow  for  the  assignment  of  physically  useful  properties  to  the 
functions  and  f2.  Specifically,  f2  is  chosen  so  that  the  second 
term  is  negligible  during  short  time  histories.  Thus,  f^  is  selected 
so  that  Eqs.  (1)  and  (7)  will  reproduce  the  family  of  stress- 


strain  curves  run  at  high  strain  rates;  the  function  £2  is 
neglected  in  this  situation.  Conversely,  both  and  £2  are 
required  to  predict  secondary  creep.  Mathematically,  £2  is  a 
perturbation  term  that  must  be  added  for  long  time  response. 

The  main  thrust  of  this  work  is  directed  toward  establishing 
representations  for  f^^  and  £2,"  and  developing  the  relationships 
between  these  representations  and  the  experimental  response  of 
the  material.  The  first  step  in  achieving  this  goal  is  to 
recognize  that  for  any  experiment,  the  stress  and  inelastic  strain 
rate  histories  can  be  determined  from  the  experimental  data. 

Then,  by  inverting  Eq.  (7)  , 


Zd-u)  =  a(t) 


2n 

n+1 


/3  eP(t) 


l/2n 


the  history  of  the  state  variable  product  Zd-o)  can  be  calculated 
for  each  particular  experiment.  This  provides  a  powerful  techniaue 
to  directly  determine  the  functions  f^,  £2  and  10. 
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SECTION  IV 


EVALUATION  OF  TENSILE  RESPONSE 

The  coefficient  n  in  Eq.  (7)  controls  the  strain  rate 
sensitivity  of  the  model  with  increases  in  n  corresponding  to  a 
decrease  in  strain  rate  sensitivity.  The  properties  of  Eq.  (7)  are 
such  that  for  each  choice  of  n,  a  family  of  rate  dependent  stress 
strain  curves  can  be  predicted  from  a  single  representation  of  Z. 

Thus,  the  family  of  Z  histories  calculated  from  Eq.  (10)  using 
experimental  data  can  be  reduced  to  a  single  response  curve  Z  with 
the  correct  value  of  n.  The  relationship  Z(W^)=  Zq  +  f (W^)  is  valid 
for  this  analysis  provided  the  time  duration  of  the  tensile  tests  are 
short  compared  to  the  creep  response  and  material  damage  is  not 
significant. 

The  response  of  Rene  95  at  650“C  published  in  Ref.  [5]  was  used 
for  this  analysis.  This  work  contains  a  family  of  four  tensile 

_  3  -4 

response  curves  run  at  strain  rates  of  8.3  x  10  ,  8.3  x  10  , 

“6  •’1 

8.3  X  10  and  1.6  x  10  sec  .  The  observed  tensile  response  at 

the  highest  strain  rate  had  a  0.2%  yield  stress  of  about  1200  MPa  and 

almost  no  strain  hardening  with  the  ultimate  stress  less  than  1300  MPa. 

The  other  tensile  response  curves  were  similar  with  a  lower  value  of 

yield  stress.  From  a  family  of  creep  curves,  it  was  observed  that  for 

the  highest  value  of  creep  stress,  1200  MPa,  the  corresponding  minimum 

creep  rate  was  2.7  x  10**^  sec^.  Thus  it  was  concluded  that  the  two 

-  3  -  4  -1 

tensile  curves  run  at  8.3  x  10  and  8.4  x  10  sec  had  a  negligible 
amount  of  creep  and  could  be  used  for  evaluating  the  rate  sensitivity 
constant,  n.  Unfortunately,  the  highest  rate  is  only  a  factor  of  ten 
greater  than  the  next  rate,  and  the  resulting  tensile  curves  crossed. 
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Thus,  we  used  the  value  of  n  =  3.2  reported  in  the  literature, 

Ref.  [7],  for  Rene  95  at  650°C. 

The  history  of  Z  as  a  function  of  plastic  work  was  established 

-4 

with  n  =  3.2  for  the  tensile  test  run  at  a  strain  rate  of  8.3  x  10 
sec  The  state  variable  history  is  shown  in  Fig.  1.  The  curve 
appears  to  be  of  an  exponential  form  with  a  saturation  value  = 

318  MPa.  The  condition  of  a  saturation  value  of  Z  agrees  with  the 
observed  behavior  that  the  stored  energy  of  cold  work  increases  with 
deformation  at  a  decreasing  rate.  A  further  investigation  showed 
that  2,n(Zj^-Z)  was  very  close  to  a  linear  function  of  the  plastic  work, 
W^.  Thus  the  representation 

Z  =  f^  +  Cq  =  Zj^  -  (Z^  -  Zq)  exp(-mwP)  (11) 

was  estciblished  from  the  data  with  the  constants  Z^  =  16  82  MPa  and 
m  =  0.37  MPa  ^  determined  using  a  least  squares  analysis.  The  initial 
hardness  Zq  was  determined  by  extrapolating  !ln(Z^-Z)  to  zero  plastic 
work.  Data  could  not  be  accurately  determined  for  small  values  of 
plastic  work;  however,  it  is  expected  that  values  of  Z  exist  below 
Zq  as  indicated  by  the  dotted  line  in  Fig.  1. 

The  first  term  of  the  evolution  equation,  Eq.  (9) ,  for  tensile 
response  can  be  obtained  by  differentiating  Eq.  (11)  to  obtain 

Z  =  m(Zj^  -  Zg)  exp(-rawP)  =  f^(wP)  W^.  (12) 

An  alternative  formulation  for  the  tensile  date  was  established 
by  Bodner,  Ref.  [7],  by  combining  Eqs.  (11)  and  (12)  to  get 

Z  =  m(Zj^  -  Z)  wP  =  f[(Z)  wP.  (13) 
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Even  though  f^(W^)  and  f^(Z)  are  equivalent  for  tensile  response 
alone,  the  two  representations  for  f^  give  two  different  forms  for 
the  evolution  equation  with  different  properties.  In  this  presen¬ 
tation  we  elect  to  investigate  the  consequences  of  using  the  second 
representation,  Eq.  (13). 

In  closing,  it  should  be  noted  that  Eq.  (11)  is  different  from 
the  representation  reported  in  Ref.  [5]  for  the  same  material.  Both 
models  essentially  the  same  value  through  the  range  of  interest; 
however  the  asymptotic  behavior  of  Eq.  (11)  appears  to  be  more 
consistent  with  the  observed  response  of  the  material. 
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STATE  VARIABLE  ,  Z  (M  Pa) 
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Figure  1.  Internal  state  variable  history  for  one  tensile  test  and  eijdu  eree]' 

tests.  The  representation  for  tensile  function  r|(w'’)  +  is  cemp.ins 
to  the  experimental  data. 


SECTION  V 

EVALUATION  OF  THE  CREEP  RESPONSE 
During  secondary  creep  (when  the  strain  rate  is  approxi¬ 
mately  constant)  the  value  of  Z  must  be  constant  if  Eq.  (7)  is 
to  be  satisfied.  Thus  for  secondary  creep,  Eq.  (9)  becomes 

Z  =  0  =  fj^'  (Z)  wP  -  f2(Z)  .  (14) 


The  stationary  value  of  Z  corresponds  to  a  stable  material 
microstructure  during  the  secondary  creep  regime. 

In  general  fj^'(Z)  is  known  from  the  tensile  data  and  the 
stationary  value  of  Z  can  be  calculated  from  Eq.  (9)  using  the 
creep  response  data.  The  Z  histories  from  eight  creep  tests 
are  also  shown  in  Fig.  1.  It  is  interesting  to  note  that  for  the 
duration  of  each  test,  Z  is  nearly  constant.  This  occurs  because 
the  amount  of  primary  creep  during  a  typical  test  is  very  small 
compared  to  amount  of  secondary  creep. 

To  start  to  establish  a  representation,  let  us  assume  the 
Bodner  representation  for  f2;  i.e. 


where  Z2  is  constant  defining  the  minimum  value  of  Z.  If 

*  *  D 

Eq.  (15)  is  valid,  then  )ln(f^  (Z)  V- )  must  be  a  linear  function 
of  Jln(  (Z  -  Z2) /Zj^)  for  Eq.  (14)  to  be  satisfied.  The 
experimental  data  was  reduced  and  plotted  as  shown  in  Fig.  2. 

The  value  Z2  =  1300  MPa  was  chosen  (arbitrarily)  so  that 
Z  -  Z2  is  always  positive. 
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Many  of  the  creep  properties  of  Rene  95  are  summarized  in 
Fig.  2.  First  notice  that  the  data  is  bilinear  with  a  sharp  change 
in  slope.  The  constants  for  each  region  were  calculated  using  a 
least  squares  analysis  and  it  was  found  that  the  curves  intersected 
at  Z  =  1682  MPa;  a  value  that  corresponds  to  a  creep  stress  of 
approximately  980  MPa.  Since  r  =  12.0  and  r  =  1.34  for  stress  above 
and  below  980  MPa,  it  appears  that  different  physical  properties 
control  deformation  in  the  "high  stress"  and  "low  stress"  domains. 
Further  observe  that  the  intersection  point,  Z  =  1682  MPa,  is  the 
initial  value  of  hardness,  Z^,  determined  in  the  tensile  tests. 

Thus  the  two  response  domains  are  segregated  by  the  apparent  initial 
value  of  hardness  Zq. 

Next  recall  that  Zq  was  estimated  by  extrapolation  to  zero 
plastic  work;  but  it  was  expected  that  values  of  Z  below  Zq  existed 
at  very  small  values  of  plastic  work.  The  creep  data  gives  further 
evidence  that  Zq  is  not  the  initial  value  of  Z.  Equation (7)  requires 
the  initial  value  of  Z  in  a  creep  test  must  be  below  the  stationary 
value  during  secondary  creep  in  order  to  predict  a  primary  creep  rate 
greater  than  the  secondary  creep  rate.  This  indicates  that  the  initial 
rate  of  hardening,  Z,  is  very  high  for  values  of  Z  below  Zq  =  1682  MFa. 
This  further  confirms  the  existence  of  different  deformation  properties 
being  active  above  and  below  a  stress  of  980  MPa. 

It  is  difficult  to  accurately  model  the  total  creep  and  tensile 
response  of  Rene  95  when  the  response  data  is  not  available  for  very 
small  values  of  plastic  work.  However,  a  reasonable  approximation  can 
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be  constructed  by  observing  that  the  low  stress  data  in  Fig,  2  is 
nearly  constant  when  compared  to  the  high  stress  data.  Let  us  set 

=  Zq  and  neglect  changes  in  Z  below  1682  MPa.  In  this  case  r  =  3.0 
and  A  =  5. 34  sec 

If  damage  is  assumed  to  be  the  major  mechanism  that  produces 
tertiary  creep,  an  estimate  of  the  value  of  the  damage,  lu,  at  failure 
can  be  obtained  from  the  creep  data.  A  total  evaluation  of  the  damage 
variable  is  not  included  since  the  creep  curves  were  established  using 
constant  load  rather  than  constant  stress  conditions.  However,  after 
adjustment  for  area  reduction,  it  was  found  that  the  value  of  Z(l-^) 
decreased  about  one  to  four  percent  between  secondary  creep  and  failure. 
It  follows  the  value  of  oj  at  failure  can  be  expected  to  be  very  small 
in  most  cases.  Thus,  the  popular  approach  of  assuming  uj  =  1  at 
failure.  Refs.  [15,  16],  is  not  applicable  to  Rene  95.  This  indicates 
that  failure  occurs  due  to  a  mechanism  other  than  the  Kachanov  damage 
parameter.  In  Ren4  95  this  probably  unstable  crack  growth. 


f 

I 
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Figure  2.  Bilinear  representation  for  the  stationary  values  of  the 
internal  state  variable  during  secondary  creep. 


SECTION  VI 

COMPARISON  OF  THE  MODEL  AND  DATA 


The  hardening  law  for  Rene  95  at  650°C  can  be  written  by 
combining  Eqs.  (9) ,  (12)  and  (14)  to  get 


7  z  -n  Z  -  Z  r 

f-  =  m  [1  -  f-1  wP  -  A  (-5—^)  (16) 

with  Zq  the  initial  value  Z.  This  can  be  used  in  combination 
with  Eqs.  (1)  and  (7)  to  calculate  the  response  of  the  material 
assuming  oi  =  0.  The  constants  developed  for  Rene  95  are: 


4  -1 

D  =  10  sec  (assumed) 
o 

n  =3.2 
Zq  =  1682  MPa 
=  2193  MPa 
E  =  1.77  X  10^  MPa 


m  =0.37  MPa 
=  1682  MPa 

A  =5.34  sec 

r  =  3.0 


where  E  is  the  elastic  modulus.  The  final  formulation  for  tVie 

hardening  rule,  Eq.  (16)  is  identical  to  that  proposed  by 

Bodner,  Ref.  [7].  The  constants  Z^,  Z^  and  m  that  were  determined 

from  the  tensile  test  are  very  close  to  those  reported  by  Bodner. 

The  constants  Z2  and  r  developed  from  the  creep  data  are 

different.  The  constant  Z2  =  Z^  was  used  as  an  approximation. 

The  above  equations  were  used  to  predict  the  tensile 
* 

response  of  Rene  95  at  constant  stress  rate.  Recall  the  constants 
were  developed  from  a  tensile  test  run  at  constant  strain  rate. 

The  predictions  are  shown  with  two  test  results  in  Fig.  3.  The 
agreement  between  theory  and  experiment  is  acceptable.  Observe 
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that  the  stress,  a  =  980  MPa,  corresponding  to  intersection  point 
in  Fig.  2,  is  very  close  to  the  yield  stress. 


The  creep  response  of  Rene  95  was  also  calculated  for  three  stress 
levels,  903  MPa,  1034  MPa  and  1206  MPa.  These  three  stresses  are  at 
the  high  and  low  ends  of  the  two  curves  in  Fig.  2  and  also  the 
intersection.  The  calculated  curves  and  experimental  data  are  shown 
in  Fig.  4.  The  minimum  creep  rate  for  the  1206  MPa  and  1034  MPa  tests 
are  very  close;  but  the  secondary  creep  rate  for  the  903  MPa  was  under 
predicted.  This  resulted  from  the  approximation  that  was  used  for  Z 
at  the  low  stress  levels. 

The  results  in  the  previous  two  sections  suggest  that  two  different 
micromechanical  mechanisms  are  active  at  high  and  low  levels  of  stress. 
This  effect  can  also  be  seen  in  the  creep  response  curves  displayed  in 
Fig.  4.  The  stress  increments  of  131  MPa  and  172  MPa  between  the  three 
creep  tests  are  approximately  equal.  However,  the  minimum  creep  rate 
changes  by  factors  of  four  and  140,  respectively.  This  difference  in 
secondary  creep  rate  reflects  the  consequences  of  a  possible  change  in 
micromechanical  mechanism. 

In  closing,  it  is  important  to  point  out  that  the  above  results 
show  how  all  phenomenological  theories  must  be  used  with  caution. 

That  is,  the  value  of  the  parameters  that  are  determined  for  a  material 
under  a  given  range  of  conditions  may  not  apply  in  another  range  of 
conditions.  For  example,  if  the  value  of  12.0  assigned  to  r  for  the 
high  stress  region,  the  predicted  response  in  the  low  stress  region 
would  be  poor.  However,  the  above  system  of  equations  could  be  modified 
to  account  for  the  physical  change  in  the  response  and  extrapolation 
in  both  the  high  and  low  stress  regions  may  be  possible.  Of  course, 
the  new  formulation  should  be  experimentally  verified  before  application. 
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Figure  3.  Calculated  and  experimental  tensile  response  of  Rene  95 
at  650°C  under  constant  stress  rate  control. 
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Figure  4.  Calculated  and  experimental  creep  response  of  Rene  95  at  650°C  at  tiigli, 


low  and  intermediate  values  of  stress. 
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